On Classical Electromagnetic Fields 

II. RAYS: THE EIKONAL TREATMENT OF GEOMETRIC OPTICS 6 

Since ancient times, the notion of ray or beam propagation has been one of the most 
enduring and fundamental concepts in optical physics. As a zeroth order approximation we 
might consider a plane wave to be a model of a beam and its propagation vector to be a 
model of a ray. This is a reasonable start, but it is a much too restricted view and we can 
do much better. What we need is a solution to Maxwell's equations which is like a plane 
wave, but limited in spatial extent. One approach, the simplest, is called variously ray, 
Gaussian or geometric optics. 

A Ma xwellian Derivation of the Eikonal Equation: 



To fully understand geometric optics in the context of Maxwell's equations 
writing the electric and magnetic fields as pseudo-simple waves — viz. 



we start by 




[Il-la] 




[Il-lb] 



where ko= ©.^ii^ =(0 /c 



It is assumed that 

e(f,(o) and h(r,co) 

are weak functions of position . 

phase function S(r,(o) is the spatially varying phase of the pseudo-simple 
the cases of pseudo-plane waves and pseudo-spherical waves the phase 
given, respectively, by 



wave. For 



function is 



The scalar 



koS(r,(o) = A:k^-l-3;k^,-l-zk, 



[ n-2a ] 
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See, for example. Max Born and Emil Wolf, Principle of Optics, Pergamon Press (1986), Chapter 3. 
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and koS(r,a)) = + 7 [ II-2b ] 

We now substitute these pseudo-simple wave expressions (i.e. Equations [ II-l]) into 
Maxwell's equations to obtain 

exp[-/koS (?,«)] I V xe (f,co)- iko VS(f,a))xe (f,a))|=-/|i.oCk(,h (f,a)) exp[-ikoS(f,a))] 

exp[-jkoS(f ,co)]|vxh(f,(o)- ikg VS(f,(o)xh(f,(o)|= i8(f ,(o)ckQe(f,co) exp[-jkoS(f ,(o 

exp[-jkoS(f ,(o)]|v |E(f,(o)e(f,(o)]- jk(,E(f,(o) VS(f,(o) e(f,(o) j= 

exp[- jkoS(f ,0))] {v •h(f,(o) - ik^ V S(f ,(o) • h(f ,«) } = 

Rearranging, we obtain 

VS(r,a))xe(r,a))-|a,oC h(r,a))= [/kg] ^Vxe(r,a)) [n-4a] 

VS(f ,0)) X h(f ,0)) +8(f,a)) ce(f ,0)) = [ik^]^ V xh(f,a)) [ II-4b ] 

VS(r,to) (r,a))e(r,a))] = [ik^y' V [8(r,a))e(r,to)] [ n-4c ] 

V S (r,a)) • h (r,a)) = [ i k^ ] V h (r,©) [ II-4d ] 

In the ray, Gaussian or geometric approximation we assume that we may neglect 
the RHS's of these equations. To get something useful we multiply through the first 
equation (i.e. Equation [ II-4a ]) as follows: 

[|ioc]"'vs(f,co)x{ Equation [ II-4a ] } [n-5a] 
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[iIqc] ^VS(r,a))x{vs(r,a))xe(r,a))-|ioCh(r,(o)}«0 [II-5b] 

Applying the "abc = bac - cab" rule^ we obtain 
[^lo c]"'{ V S(f,a)) [ V S(f,a)) • i(f,to)] 

-e(f,(o) |vs(r,co)|'-|iocVS(f,a))xh(f,(o)[ ==0 



[n-5c] 



which becomes upon substitution from the second Equation [ II-4b ] 
[|ioc]"'{vs(f,(o)[vs(f,(o)e(f,(o)] 

- ,0)) |v S(r,o))| i+8(r,o)) ce(f ,co) ~ 



[ II-5d ] 



From Equation [ II-4c ], we see that the first term vanishes in the geometric 
approximation — i.e., if we neglect the term [^k^ ] V je (r,co) e(r,(o)] . Therefore, 
for non- vanishing e(f,03) we obtain the following reduction of Maxwell's 
equations : 

VS(f,(jo)|' ==E(f,(o)|io c' = n'(f,co) [ II-6 ] 

where '^(f ,co) is the index of refraction. More explicitly, we may write an equation for 

a "ray vector" ~ i.e. the tangent to a space curve orthogonal to the surfaces of constant 
S(r,(o) 



Again using ax(b X c) = b(a • c) - c(a • b) . 
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VS(r,a)) = n(f ,0)) s = ^(r,©) — 

ds 



We illustrate the geometric relationships below: 



[11-7] 





Reference poin't 



grad S(r^w) 



We may now derive the all important eikonal equation. To that end, we first take a derivative 
along the ray direction - viz. 



— [vs(f,co)l= — n(r,(o) — 



ds 



ds 



[ II-8a ] 



However, from the definition of the grad operator we know that 

i/[vs(?,a))] = dY-V [vs(?,a))] 



so that 



d 



df 



f • vfvs(?,co)l = 1-^ VS(f,a))|. Vfvs(f ,co)l = -f n(f ,«) ^ 
ds ^ ^ [«{r,toj J ^ dsl ds 



[ II-8b ] 
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or 




1 



v{vs(f,a))VS(f,a))} 




1 



v{n'(r,a))} = Vn(r,a)) 



J r \ dr 
= — n r,a) — 
dsl ^ ^ ds 



[ n-8c ] 



Thus we have obtain the eikonal ^ equation for the ray vector ~ viz. 



d \ dr 

— n[r,G)) — = 

ds\_ ds _ 




[11-9] 



First Application of the Eikonal Equation: Mirages 

Air adjacent to a hot surface rises in temperature and becomes less dense. Thus over a 
flat hot surface, such as a desert expanse or a sun drenched roadway, air density 
locally increases with height and the average refractive index may be approximated 
by a simple linear variation of the form 



where x is the vertical height above the planar surface, ng is the refractive index at 
ground level, and K is a positive constant. 

We may use the eikonal equation to find an equation for the approximate ray 
trajectory — i.e. an equation for ray height x as a function of ground distance z -- of a 
hght ray launched from a height Xq and at an angle with respect to the surface of the 

earth. 



n {x)= {1+ K^} 



[11-10] 



The eikonal (from the Greek: eiKOSV means image) was introduced in 1895 by H. Brans. 
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Therefore, 



d_ 
ds 



/- ^ dr 
n r,0) — 

' ds. 



= V/i(f,a)) 



d^x 1 



dz n[x,Gi) dx 



d , . 
— nyxisy) 



[n-lla] 



or from Equation [ II- 10 ] 



d^x 
Hz' 



[Il-Ub] 



Thus, the ray trajectory is given by 



f(z) = j^-^ z^ + tanSg z+Xqj x+z 



[11-12] 



Ray trajectories diverted by a hot surface 




20 40 60 SO 100 

Distance along hot pUnar surface 
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Second Application of the Eikonal Equation: The "ABCD" Ray 
Matrices - A Systems Approach to Optics 



1. Uniform dielectric medium — i.e. Mj,(xi) is a constant so that — 

as 



dr 

. ds . 



0. 



Thus, the ray must be a straight line which may be written r = 5 a + b . 

In the two-dimensional paraxial approximation, we assume that s~ z and write 



J dp 

Pout =Pm + L — 

dz 



and ^ 

dz 



dp 

dz 



[11-13] 



where p = p0=xx + yy 



Pin 




^P'out 

Pout 


_- 


^ in 

^in 



1 


— 


^ out 



We may write results of this sort in the form of the famous and highly useful ray 
transform or ABCD matrix -- viz. 



Pout 




A 


B 




Pout. 




c 


D 





[11-14] 
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In the case of a uniforai dielectric 



Pout 




1 


L 




Pout. 







1 





[11-15] 



so that A = 1, B = L, C = 0, and D = 1 



2. A dielectric discontinuity : Starting with Equation [ II-7 ] and noting, once again, 
that curl grad { }=VxV{ } = Owesee that 

VxVS(f,a))= Vx|/i(f,a)) s}= [11-16] 

which is identical to the saltus condition on the electric and magnetic fields at a 
dielectric interface! Hence | n\r,(iij s [ continuous across the dielectric 

J tangent 

boundary so that sin6j = sinGj — i.e. Snell's law! This result in the paraxial 
approximation ( i.e., sinG » tanG = G ) may be written in ray matrix form as 



Pout" 




"A 


b' 




"pin" 




"1 




"P,n" 


Pout. 




c 


D 




.pL 




.0 "in/ "out. 




.pL. 



[11-17] 
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3. A " Thin" lenses : In passing we note that the ray matrix of a thin lens is given by 
or, perhaps more accurately, a thin lens is essentially defined by 



[11-18] 



Pout" 




'A 


B' 




"P,n" 




1 




Pin' 






c 


D 















4. Axially symmetric GRIN media : Consider the use of GRaded INdex technology to 
obtain an axially symmetric variation in the index of refraction of the form ^ 



n{p)= 



r 



l-A 



rp 



'1 



[11-19] 



A note on GRIN technology : In GRIN technology one builds up a glass rod with a specific radial index of 
refraction distribution by fusing a sequence of coaxially arranged glass tubes with appropriate index and diameter 
as illustrated in the following: 

Coaxial dielectric (glass) tubing 
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Within such a GRIN rod, we write f = pp + z z for the ray coordinates and 

V«(f) = (3 — n(p) for the index variation. Using the eikonal equation — i.e. 
dp 

Equation [ II-9 ] ~ in the paraxial approximation, we find 



dz 



[ II-20a ] 



or 



d . s J r . . dp~\ 

nip) ~ — I n p — I 

dp dz[_ ^ ' dz\ 



[ II-20b ] 



Therefore 



— 2 P ' / \ 
dz n\p) dp 



^ :^"(p)=-^i4«(p)] 



[ n-20c ] 



or 



:p- 



dz^ n{p) dp y 



1- A 



fpTll 1 r mAfpf-'l 
\a J \ n\p) a \a J 



niA 



a 



/ \m— 1 
P 



[ II-20d ] 



Doubtless, the simplest and most valuable instance is m = 2 ~ i.e. what is usually 
called parabolic or quadratic material ~ wherein 



d^ 2A 

— P- 

a 



dz 



^pf 2A 

- =- — P=-K P 

Va y a 



so that 



P (z) = Pin cos(kz) + ^ sin(Kz) 



[11-21] 
[ II-22a ] 
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P\z) =-PinK sin(Kz)+p; cos(kz) 
In tenns of a ray transform matrix 



Pout 




A 


B 




J^out. 




c 


D 





r cos(kz) k ^sin(Kz)1 
[-Ksin(Kz) cos(kz) J 



where K= yjlAfi 



Ray trajectories confined in a GRIN rod. 



[ II-22b ] 



[ 11-23 ] 



Imaging Distance 
K AZ = TL 




Focusing Distance 
K AZ = n/2 
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Alternative (Hamiltonian) Derivation of Eikonal Equation: 

FERMAr S PRINCIPLE 

Like most laws of physics, the equations of geometric optics can be derived from a 

variation principle. In this context the variation principle is called the Fermat 
principle which states that a ray always chooses a trajectory that minimizes the optical 
path length - viz. 

J n {x,y, z) ds = minimum [ 11-24 ] 

Pi 



See, for example, Dietrich Marcuse, Light Transmission Optics, Van Nostrand Reinhold (1972). 

11 More precisely, the path must be a local extremum and in rare cases may, in fact, be a maximum. See R. Y. 
Luneberg, Mathematical Theory of Optics, University of California Press, Berkeley and Los Angeles (1964). 
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where the line element, ds , is measured along a ray and the two end-points Pi and P2 
are fixed in space. 12 Analysis of the variation problem is facilitated by choosing the 
projected coordinate z as the new variable of integration. Accordingly, 

ds= -Jdx^ + dy^ + dz ' = yjl + x^ + y^ dz , [ 11-25 ] 

, dx , dy 

where x = — and y = —, Fermat's variation principle is transformed into the more 
dz dz 

familiar Lagrangian form ~ viz. 

Pi 

I L{x, y, x', y')dz = minimum [ 11-26 ] 

Pi 

where L{x,y, x',y') = r^x,y,z) ■^l^^x^^'^^Y^ . [11-27] 

The minimization procedure is then well-known in the variational calculus and 
leads to the famous Euler-Lagrangian equations ~ i.e. 

d dL dL „ 

— —,- — =0 [n-28a] 

dz ox ox 

d dL dL ^ 

— —,-—=0 [II-28b] 
dz oy ay 

When applied to the Fermat Lagrangian, as defined in Equation [ 11-27 ], these 
equations yield 



1^ From Equation [ 11-7 ] we see that 

S{F,)- S{F,)= J n{x,y,z)ds . 
Pi 
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= ^\ + x^ + y^— [II-29a] 



ny' r Ti^n 

=d===^l + x +y — . [II-29b] 



dz ^i+x'^+y^ ^ dy 

Using Equation [ 11-25 ] we see that the Euler-Lagrangian equations may be 
expressed in the vector form as 

[jT dx 1 d I dy 1] _ \ d n d n] nSO 

which is precisely the content of Equation [ II-9 ] -- QED. 

Hamiltonian Formulation of Ray Optics 

The analogy between ray optics and particle mechanics is most striking when the 
equations of ray optics are expressed in Hamiltonian form.i^ To that end, we define the 
generalized momentum which is canonically conjugate to x and y by the vector 
equation 

The Hamiltonian is then define in terms of the generalized momentum by the relation 

H(x,y,p,,Py)= p,x+Pyy-L{x,y,x,y) . [11-32] 

With the assumed functional dependence of the Hamiltonian, we form the derivatives 



1^ The formal theory of optical systems was developed by Sir W. R. Hamilton in 1828-37. 
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dH , dx' dy' dL dx' dL dy' 

^=^+P.^ + Py^ T~^T~ [II-33a] 

op, dp^ dp^ dx dp^ dy dp^ 

dH dx' , dy' dL dx' dL dy' ^ „ , 

^=P^^+y+Py^-^'^ T-^T^ • [II-33b] 

dPy dpy dpy dx dpy dy dp^ 

Given the definitional relationships embodied in Equation [ 11-31 ] we see that these 
expression reduce to one set of Hamilton's equation ~ viz. 



\dx_ dy\_ \dH_ dH [ 
[dz'dzj [dp/ dp J' 

The other set of Hamilton's equation ~ viz. 

\ dp^ ^1 

{dz'dz] [ 3x ' dy 



[11-34] 



[ 11-35 ] 



follow directly from the Euler-Lagrangian equations ~ i.e. Equations [ II-28a ] and 
[ II-28b ] ~ and the definitions embodied in Equation [ 11-31 ]. Using the Fermat 
Lagrangian we see that 

r 1 f 3L 3L 1 [ nx' ny' 1 

{p.,Pyh\ = / . , [11-36] 

[dx dy ] yj\+x'^+y^ + + j 
and consequently that we may solve for { x',y' } in terms of { p,,Py } as 

^ P P ^ 

{-VIH ■ ^ , \ [11-37] 

[ -P.-Py V" -P.-Py J 
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Substituting into Equation [ 11-32 ], we find an expression for the Fermat or ray 
optics Hamiltonian ~ viz. 

which resembles the mechanical Hamilton of a relativistic particle ~ i.e., 



But the analogy is even stronger in the paraxial approximation where the Hamiltonian is 
approximated by an expression which is identical in form with the Hamiltonian of a 
non-relativistic particle ~ viz. 

rr I Px+Py Px+Py r XT 



when and « {n) 



Applying the quantization rules of quantum mechanics to these Hamiltonians, we can go full circle and recover 
wave optics from ray optics. Equation [ 11-38 ] leads directly to the equivalent of the relativistic Klein-Gordon 
equation while the equivalent of the nonrelativistic Schrodinger equation follows directly from Equation [ 11-39 |. 
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